Gaps between jets at hadron colliders in the next-to-leading BFKL framework 
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We investigate diffractive events in hadron-hadron collisions, in which two jets are produced and 
separated by a large rapidity gap. In perturbative QCD, the hard color-singlet object exchanged in 
the t— channel, and responsible for the rapidity gap, is the Balitsky-Fadin-Kuraev-Lipatov (BFKL) 
Pomeron. We perform a phenomenological study including the corrections due to next-to-leading 
logarithms (NLL). Using a renormalisation-group improved NLL kernel, we show that the BFKL 
predictions are in good agreement with the Tevatron data, and present predictions which could be 
tested at the LHC. 
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The discovery of difTractive interactions in hard processes at hadron colliders has triggered a large effort devoted to 
understand rapidity gaps in processes with a hard scale. Since events with a gap between two jets were observed in pp 
collisions at the Tevatron about 10 years ago many theoretical works have investigated these processes in QCD. 
\Q • While describing diffractive processes in QCD has been a challenge for many years, one should be able to understand 
£N| , hard diffractive events with perturbative methods. 

In a hadron-hadron collision, a jet-gap-jet event features a large rapidity gap with a high— St jet on each side 
(Et^A-qcd)- Across the gap, the object exchanged in the t— channel is color singlet and carries a large momentum 
O ■ transfer, and when the rapidity gap is sufficiently large the natural candidate in perturbative QCD is the Balitsky- 
Fadin-Kuraev-Lipatov (BFKL) Pomeron 0. Of course the total energy of the collision y/s should be big (y/s ^> Et) 
in order to get jets and a large rapidity gap. 

To compute the jet-gap-jet process in the BFKL framework, the problem of coupling the BFKL Pomeron to partons, 
i— 1 1 instead of colorless particles in the standard approach, has to be solved first. Indeed, one usually uses the fact that 
impact factors vanish when hooked to gluons with no transverse momentum, a property of colorless impact factors. 
For instance, this is what allows to turn the Feynman diagram calculation into a conformal invariant Green function 
Q. However, with colored particules, this BFKL Green function cannot be used and should be modified accordingly. 
' The Mueller- Tang (MT) prescription [4j is widely used in the literature. 
. On the phenomenological side, the original MT calculation is not sufficient to describe the data. A first try to 
improve it was attempted jf| by also taking into account the soft interactions which in hadron-hadron collisions can 
\ destroy the rapidity gap. An agreement was only obtained if the leading-logarithmic (LL) BFKL calculation was done 
with a fixed value of the coupling constant a s , which is not satisfactory. Moreover, it is also known that next-to- 
leading logarithmic (NLL) BFKL corrections can be large. In Q, it was shown that a good description of the data 
could be obtained when some NLL corrections were numerically taken into account in an effective way [Tj , but the full 
NLL-BFKL kernel [§[ could still not be implemented. As a result, these tests on the relevance of the BFKL dynamics 
were not conclusive. 

On the theoretical side, it is known that NLL corrections to the LL-BFKL predictions can be large due to the 
appearance of spurious singularities in contradiction with renormalization-group requirements. However it has been 
realised 0, [13] that a renormalisation-grou p im proved NLL regularisation can solve the singularity problem and lead 
to reasonable NLL-BFKL kernels (see also [11] for different approaches). This motivates our phenomenological study 
of jet-gap-jet events in the NLL-BFKL framework. One goal is to motivate further measurements at the Tevatron 
and the LHC. 
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Such NLL-BFKL phenomenological investigations have been devoted to the proton structure function [12j , forward 
jet production in deep inelastic scattering [ill, [lij], and Mueller- Navelet jets in hadron-hadron collisions [HI, EE]- 
While for the structure function analysis the NLL corrections did not really improve the BFKL description, it was 
definitively the case in the forward jet analysis, and they play a non-negligible role for the Mueller-Navelet jet 
predictions. Therefore, it is natural to investigate what happens for jet-gap-jet events, and what is the magnitude of 
the NLL corrections with respect to the LL-BFKL results for this observable. 

The plan of the paper is the following. In section II, we introduce the phenomenological NLL-BFKL formulation 
of the jet-gap-jet cross section. In section III, we compare the LL- and NLL-BFKL calculations, and confront them 
with the Tevatron measurements. In section V, we present predictions for the jet-gap-jet cross section at the LHC. 
Section V is devoted to conclusions and outlook. 



II. THE JET-GAP-JET CROSS SECTION IN THE BFKL FRAMEWORK 



In a hadron-hadron collision, the production of two jets with a gap in rapidity between them is represented in Fig.l, 
with the different kinematic variables. We denote y/s the total energy of the collision, px and — pt the transverse 
momenta of the two jets and x\ and x-i their longitudinal fraction of momentum with respect to the incident hadrons 
as indicated on the figure. The rapidity gap between the two jets is Ar] = \n.{xiX2s/p T ). We write the cross section in 
the following form 

d(T PP^XJJY d<7 99 ^ 99 

d Xl d X2 dpl = S f«ffo>PT)f«fi*2,Pr)-^, (D 

where the functions feff(x,p T ) are effective parton distributions that resum the leading logarithms log(p^/AQ CD ). 
They have the form 

C 2 

f eff {x,^ 2 ) = g{x,^ 2 ) + — |(g(a;,// 2 ) + q(x,fi 2 )) , (2) 

c 

where g (respectively q, q) is the gluon (respectively quark, antiquark) distribution function in the incoming hadrons. 
Even though the process we consider involves moderate values of x\ and x 2 and the perturbative scale p\ 3> Aq CD , 
which we have chosen as the factorization scale, the cross section fl} does not obey collinear factorization. This is due 
to possible secondary soft interactions between the colliding hadrons which can fill the rapidity gap. Therefore, in JTJ), 
the collinear factorization of the parton distributions / e / / is corrected with the so-called gap survival probability S : 
since the soft interactions happen on much longer time scales, they are factorized from the hard part da"~^ 99 / ' dp\. 
This cross section is given by 

da 99 ^ 99 1 , „ /A 2 ,,2 

— = — \A(A V , P 2 T )\ (3) 
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in terms of the gg — ► gg scattering amplitude A(Ar],p 2 ^). The two measured jets are initiated by the final-state gluons 
and we are neglecting possible hadronization effects. 

In the following, we consider the high energy limit in which the rapidity gap Arj is assumed to be very large. The 
BFKL framework allows to compute the gg — > gg amplitude in this regime, and the result is known up to NLL 
accuracy. Since in this calculation the BFKL Pomeron is coupled to quarks or gluons, the BFKL Green function 
cannot be used as it is and should be modified. The transformation proposed in [3] is based on the fact that one 
should recover the analiticity of the Feynman diagrams. It was later argued that their prescripion corresponds to a 
deformed representation of the BFKL kernel that indeed could be coupled to colored particules and for which the 
bootstrap relation is fulfilled |l8|. Applying the MT prescription at NLL leads to 

2 leN&at* ^ f d 1 [p 2 - ( 7 -l/2) 2 ]exp{q(p|) Xe// [2p, 7 , afr^A??} 
1 V,PT> C F p\ J^J 2m [( 7 -1/2)2 _(p_ 1/2)2] [(7 _ 1/2)2 i m U 

with the complex integral running along the imaginary axis from 1/2— ioo to 1/2+ioo, and with only even conformal 
spins contributing to the sum [l9(. The running coupling is given by 



a(p 2 T )=a s (p 2 T Wir=[blo S ( P yAi 1CD )]- 1 , ft = H^L_?*> . (5 ) 




FIG. 1: Production of two jets separated by a large rapidity gap in a hadron-hadron collision. The kinematic variables of 
the problem are displayed, s is the total energy squared of the collision, pr (ni) and — pr (772) are the transverse momenta 
(rapidities) of the jets and xi and X2 are their longitudinal momentum fraction with respect to the incident hadrons. Ajy is the 
size of rapidity gap between the jets. 



Let us give some more details on formula (fj| . The NLL-BFKL effects are phenomenologically taken into account by 
the effective kernels Xeff{p, 7, <5). For p — 0, the scheme-dependent NLL-BFKL kernels provided by the regularisation 
procedure Xnll(j,u) depend on u>, the Mellin variable conjugate to exp(A?7). In each case, the NLL kernels obey a 
consistency condition [|[ which allows to reformulate the problem in terms of Xeffd, &)• The effective kernel Xeff{~f, Oi) 
is obtained from the NLL kernel xnll{"1,u) by solving the implicit equation Xeff = Xnll{i,& Xeff) as a solution 
of the consistency condition. 

In [IE, mi, the regularisation procedure has been extended to non-zero conformal spins and the kernel xnll(p, 7, 
was obtained, the formulae needed to compute it can be found in the appendix of [15j (in the present study we shall 
use the S4 scheme in which xnll is supplemented by an explicit a dependence, the results in the case of the S3 
scheme are similar). Then the effective kernels Xeff(p, 7, <5) are obtained from the NLL kernel by solving the implicit 
equation: 

Xeff = Xnll(p, 1, a Xeff) ■ (6) 

It is important to note that in formula (fj|), we used the leading-order non- forward quark and gluon impact factors. 
We point out that the next-to-leading-order impact factors are known [20J, and that in principle a full NLL analysis 
is feasible, but this goes beyond the scope of our study. 

By comparison, the LL-BFKL formula is formally the same as (UJ) with the substitutions 

Xe//(P,7.«) -+XLL{p,l) = 2 V ; ( 1 ) -i>( 1 -7 + -7T J -^(7+ -TT ) . a(fc 2 ) -> a = const, parameter, (7) 
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where ^{l) — d log IX7) / dr/ is the logarithmic derivative of the Gamma function. In the NLL-BFKL formula, the value 
of a is imposed by the scale p\ and the variation of as is given in Formula (5). We shall later test the sensitivity of 
our results when using Ap^ and varying A between 0.5 and 2. This is done using formula (4) with the appropriate 
substitutions (l3l. Ii~4j| 

a{p 2 T ) -> a(Xp 2 T ) + ba 2 (p 2 T )HX) , 

An -> A77 - ln(A) . (8) 



By contrast, in the LL-BFKL case that we consider for comparisons, a is a priori a parameter. We choose to fix it 
to the value 0.16 obtained in [l3| by fitting the forward jet data from HERA. This unphysically small value of the 
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FIG. 2: Ratio between the conformal-spin components p = 0, 1, 2, 3 and 4, and the total BFKL cross section at LL and 
NLL accuracy. Left plot: as a function of Et for Ar/ > 4; center plot: as a function A77 for the second leading jet with 
15 < Et < 25 GeV; right plot: as a function Arj for the second leading jet with Et > 30 GeV. 



coupling (0.16) is indicative of the slower energy dependence of the forward jet data compared to the LL-BFKL cross 
section. And in fact, the value a = 0.16 mimics the slower energy dependence of NLL-BFKL cross section, which in 
the forward jet case is consistent with the data. Therefore in both cases one deals with one-parameter formulae: the 
absolute normalization is not under control. In the NLL case, this is due to the fact the we do not use NLO impact 
factors, although the correction should be of order one. 



III. COMPARAISON WITH TEVATRON DATA 



The DO collaboration has performed a measurement of the jet-gap-jet event rate as a function of the second leading 
jet transverse energy Et = \pt\, and also as a function of the rapidity difference A77 between the two leading jets 
They requested two jets reconstructed in the DO calorimeter with Et > 15 GeV for the second leading jet 
and the following cuts on the jet rapidities: 1.9 < | r/x ,2 1 < 4.1 and 771772 < 0. The difference in rapidity between 
both jets A77 was imposed to be bigger than 4. The DO collaboration measured the ratio of that production cross 
section with a gap of at least 4 units of rapidity between the jets, to the dijet inclusive cross section. The first 
measurement performed presented the jet-gap-jet ratio as a function of the second leading jet Et, and the second and 
third measurements displayed the ratio as a function of Arj for the low— Et and high— Et jet samples respectively 
(low Et means 15 < Et < 25 GeV and high Et means Et > 30 GeV where Et means the transverse energy of the 
leading and the second leading jets). 

In the following, instead of using X\ and x 2 we will work with the rapidities of the two jets r\\ — \vl(x\^[s j Et) and 
r}i = — ln(x2y/s I Et) , and with the average rapidity y = (771 + ?72)/2 : 

E T E t da pp ^ XJJY da pp ^ XJJY 

Xl = -= exp(Ar;/2 + y) , x 2 = -= exp(A7?/2 - y) , = Xl x 2 2 ■ (9) 

^/ b V (lyUjLxijQilL/rp u,X\LLX2U'Pt 

For comparison with the different measurements, we integrate the NLL-BFKL cross section over y and A77, or y and 
Et, taking into account the cuts detailed above. 



A. Effect of non-zero conformal spins and NLL/LL comparison 



Before comparing our BFKL calculations with the Tevatron measurements, it is worth discussing the effect of taking 
into account all the conformal spins at LL and NLL which was not performed up to now. In Fig. [H we give the ratio 
of the different conformal spin components for p = 0, 1, 2, 3 and 4 with respect to the sum of all components. First 
we see that the p — component dominates the cross section, however the contribution of the p = 1 component is 
large at high Et and low A77, where it can reach up to 30%. Its Et dependence is also not negligible and varies 
between 10 to 30% for a tranverse energy between 15 and 60 GeV. It is also worth to notice that considering only the 
conformal spin p = component instead of computing the sum of all components leads in general to larger differences 
at LL than at NLL, especially at high Et- 
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FIG. 3: Left plots: LL- and NLL-BFKL cross sections as a function of the momentum transfer Et (upper plot) and the rapidity 
gap Ar/ (lower plots, 15 < Et < 25 GeV and Et > 30 GeV); we also give the cross sections when only the conformal spin 
component p — is included. Right plots: ratio of the NLL- and LL-BFKL calculations. Since both normalizations have been 
adjusted to reproduce the data, the ratios are close to 1. 



Let us now evaluate the relative NLL corrections with respect to the LL-BFKL cross section. In Fig. [31 both 
calculations are compared with their normalizations adjusted to describe the data (this will be discussed next), 
therefore the cross sections displayed on the left plots are similar. Considering the shapes of the curves, the differences 
between the NLL and LL calculations arc in general small. This is due to the small effective value used for a in the 
LL calculation. However one can see on the right plots, in which the NLL/LL ratios are displayed, that the differences 
are sizeable in the case of the Ar] dependence for high-£r jets (Et > 30 GeV). 



B. Comparison with the DO measurement 



In Fig. 4, we compare our predictions with the DO measurements. The jet-gap-jet cross section is computed using the 
LL- or NLL-BFKL formalism while the total cross section for dijets separated by an rapidity interval A77 is calculated 
using the NLOJct++ program [21] . To obtain the NLOQCD predictions of the dijet cross section, by which the 
jet-gap-jet cross section is divided, the parton-parton hard cross section is computed at next-to-leading order with 
respect to a s , and the leading and next-leading logarithms a™ log" Ej, and a™ log" -1 Ei, are resummed using the 
DGLAP equations [H| which govern the evolution of the parton distribution functions (PDFs) of the proton. We 
used the CTEQ6.1M parton distribution functions [23| . and the renormalization and factorization scale fi r and fif 
were chosen as fi r = /i/ = Eti, the transverse energy of the leading jet. 

It is worth pointing out that these calculations are performed at parton level, and the hadronization effects cannot 
be taken into account at this stage. However, the fact that the DO collaboration measured ratios of cross sections 
reduces the influence of hadronization corrections and of the choice for the jet algorithm, since most of these effects 
are the same for jet gap jet or dijet events. This is however not the case for underlying events corrections which do 
not play any role for jet-gap-jet events since a gap is observed. These corrections depend only smoothly on Et and 
Ar] in the limited range of Et of the DO measurement, and their largest effect is taken into account in our approach 
with the gap-survival probability of S = 0.1 at the Tevatron. 

Since the normalizations of our BFKL calculations are not under control, we fit them to the DO measurement for 
the jet-gap-jet cross section ratio as a function of Et- We take the BFKL LL or BFKL NLL formalism for the jet gap 
jet cross section and the NLO QCD calculation for the inclusive cross section. For the NLL-BFKL (resp. LL-BFKL) 
calculation, the normalisation is 1.0 (resp. 0.84) with a x 2 /dof of 0.9/6 (resp. 1.7/6) when one considers only the 
points with Et > 25 GeV in the fit. Interestingly enough, the normalizations are both close to 1, and our NLL-BFKL 
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FIG. 4: Scale dependence uncertainty of the BFKL NLL calculations. The scale uncertainty is evaluated by modifying the E% 
scale used by default to or 2Ej<- The effect of the scale uncertainty is of the order of 10-15%. 



prediction with leading-order impact factors has the correct normalization and does not need to be adjusted after 
all. It is worth also noticing that the normalisations for the BFKL LL or the BFKL NLL calculations are similar. 
It is quite contrary for the QCD calculation of the jet cross section when both jets are separated with a interval in 
rapidity larger than Ary. The NLO QCD result is about 20 times larger than the LO one and it would be useful to 
know the QCD cross section at NNLO to check if the higher order corrections are still large. The resulting jet-gap-jet 
cross section ratios are displayed in the upper plot of Fig. \5\ We notice that the Et dependence of the cross section 
is well described, with the exception of the lowest Et point where the measurement is lower than our calculation 
(this is why we did not include it in the fit to compute the normalisation). This lowest Et point is definitely more 
sensitive to hadronization or underlying event corrections, and this might be the reason why it is poorly described. 
We use the same normalisations to predict the rapidity dependence of the jet-gap-jet cross section ratio for the low- 
and high— Et samples and the result is shown on the lower plots of Fig. [3 We also notice a good agreement between 
our calculations and the DO measurement, without any further adjustment of the normalizations. This is a powerful 
result, the NLL-BFKL prediction has the correct normalization for the three measurements. 

In Fig. [U we give the scale dependence uncertainty of the BFKL NLL calculations. As we mention in Section II, 
the scale uncertainty is evaluated by modifying the E T scale used by default to E T /2 or 2E T . The uncertainty of the 
BFKL NLL prediction is of the order of 10-15%. 

The fact that the LL-BFKL cross section is in good agreement with the data when a fixed value of the coupling 
is considered is consistent with the results of The fit performed on the data as a function of Et showed that 
including NLL corrections in the BFKL framework improves the description of the data. This is also consistent 
with the findings of @, where some NLL corrections were effectively included in a numerical way. In this study we 
confronted for the first time predictions obtained with the full analytic expression of the NLL-BFKL kernel to the 
Tevatron data, and the result is satisfactory. 
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FIG. 5: Comparison between the DO measurement of the jet-gap-jet event ratio with the NLL- and LL-BFKL calculations. The 
NLL calculation is in good agreement with the data, without the need to adjust the normalization. The LL calculation gives 
also a good description of the data after its normalization has been adjusted, although the fit shows that the NLL description 
is better. 



IV. PREDICTIONS FOR THE LHC 



Using the normalizations obtained from the fits to the DO data, we are able to predict the cross sections at the 
LHC. In doing so, we also have to change the gap-survival probability from 0.1, the value used for the Tevatron, to 
0.03, the proper value for the LHC. The values of the NLL-BFKL cross sections are shown in Fig. [5] as a function 
of Et for different Arj ranges, and in Fig. [7] as a function of Arj for different Et ranges. In both figures, the cross 
sections are displayed on the left plots, and on the right plots they are divided by the NLOJet++ prediction for the 
inclusive dijet cross section. In Fig. [H we display the scale uncertainties of the BFKL NLL calculation by modifying 
the default E T scale to E T /2 or 2E T as a function of the Arj between the two jets (the uncertainty is the same as a 
function of jet Et)- The uncertainty of the BFKL NLL prediction is of the order of 10-20%. 

Comparing the relative contributions from the different conformal spin components, we could see that the p = 1 
and p — 2 components are large at small A77 and at large Et, and they cannot be neglected as it is shown in Fig. [9] 
At large Arj and low Et, the p = component is by far the largest one. Also, it can be noticed that the cross 
section is high enough to be measured at the LHC in special runs, collecting integrated luminosities of a few 100 
pb _1 . Nevertheless, a good energy calibration will be needed to deal with not so energetic jets which are in addition 
quite forward. Finally, a clear prediction is that the percentage of the jet-gap-jet events is higher when the interval 
in rapidity between the two jets A77 is large. 



V. CONCLUSIONS 



Within the BFKL framework at NLL accuracy, we have investigated diffractive events in hadron-hadron collisions 
in which two high— Et jets are produced and separated by a large rapidity gap Arj. Using renormalization-group 




FIG. 6: Left plot: NLL-BFKL jet-gap-jet cross section as a function of Et for three different regions in Ar] at the LHC. Right 
plot: the left-plot curves are divided by the NLO QCD dijet inclusive cross section. 



improved NLL kernels in the S4 scheme, the NLL-BFKL effects were taken into account through the effective kernel 
obtained from the implicit equation (J6]). We implemented the MT prescription to couple the BFKL Pomcron to 
colored impact factors, this allowed our phenomenological study of NLL-BFKL effects in jet-gap-jet events. We point 
out that only used leading-order impact factors, the implementation of the NLO impact factors goes beyond the scope 
of our phenomenological analysis. 

Dividing our jet-gap-jet BFKL calculations by the NLO QCD preditions for the inclusive dijet cross section, we 
showed that the BFKL predictions were in good agreement with the Tevatron data for both the Et and A?y dependence 
of the jet-gap-jet cross section ratios. In the case of the NLL calculation, adjusting the normalization was not needed 
for any of the three measurements. The LL-BFKL predictions can also describe the data with the fixed value of 
the coupling a — 0.16, and a normalization factor of order one. Still the better description was obtained with the 
NLL formulation. We presented predictions which could be tested at the LHC, for the same jet-gap-jet event ratio 
measured at the Tevatron, but for larger rapidity gaps which can be obtained at the LHC. This should provide a 
strong test of the BFKL regime. 
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FIG. 7: Left plot: NLL-BFKL jet-gap-jet cross section as a function of the rapidity gap Ar/ for two different Et regions at the 
LHC. Right plot: the left-plot curves are divided by the NLO QCD dijet inclusive cross section. 
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scale used by default to or 2E^. The effect of the scale uncertainty is of the order of 10-15%. 
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FIG. 9: Ratio between the conformal-spin components p = 0, 1, 2, 3 and 4, and the total BFKL cross section at NLL accuracy. 
Left plot: as a function of Et in different regions of A77; right plot: as a function Ar/ for a leading jet with 20 < Et < 200GeV 
and 50 < E T < 200GeV. 
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